Quaternary-singlet State of Spin-1 Bosons in Optical Lattice 
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We present the quantum ground state properties of ^^Na spinor condensates, which is confined 
in a periodic or double-well potential and subject to a magnetic dipole-dipole interaction between 
nearby wells. A novel singlet state arise in the system and can be discussed in explicit form. Caused 
by the competition between the intra-site spin exchange interactions and the inter-site dipole-dipole 
interactions, this quaternary-singlet state is a entangled state formed by at lest four particles and 
vanish the total spin. This is distinct from the direct product of the two conventional singlet pairs. 
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I. INTRODUCTION 



The Heisenberg model of spin-spin interactions de- 
fined hyH = ^X)<i>j> • is often considered as the 
starting point for understanding many complex magnetic 
structures in solids such as ferromagnetism and antifer- 
romagnetism at temperatures below the Curie tempera- 
ture. This Hamiltonian arises from the direct Coulomb 
interaction among electrons and the Pauli exclusion prin- 
ciple, with Fi the spin operator for the ith electron. 

The study of bosonic spin-spin interactions rise since 
the success in trapping a ^"^Na condensate in an optical 
potential [1], where spin degrees of freedom are liber- 
ated and it give rise to a rich variety phenomena such 
as spin domain formations , spin mixing dynamics Q , 
topological defects and so on. The properties of such 
a three-component spinor condensate were first stud- 
ied with atomic spin coupling interaction takes the form 
V{v) = (co + C2F1 • F2)5{v) by Ho[| and Ohmi Q, and 
had been implemented experimentally Q , where two dif- 
ferent spin-dependent phases exist: the so-called antifer- 
romagnetic and ferromagnetic states for ^'^Na and ®^Rb 
atomic condensates respectively. 

This spinor BECs can be confined in optical lattices, 
which offers a unique opportunity to study magnetic 
properties of matter with tunable parameters. The quan- 
tum phase transition from the superfluid phase to the 
Mott insulating state is well described by the spinless 
Bose-Hubbard Hamiltonian@, 0] and was demonstrated 
in experiments 8, 9]. With the atomic spin coupling inter- 
action involved, the Bose-Hubbard Hamiltonian is stud- 
ied in the Mott insulating regime, where phase coher- 
ence or superfluidity is lost, and atoms are localized 
with number fluctuations suppressed. Such a insulat- 
ing state represents a correlated many-body state of 
bosons and the calculations of the magnetic properties 
within spin-exchange interactions have been carried out 
sufficiently. 10- 16] . 

In addition to the spin-exchange interaction, there ex- 
ists another important type of magnetic interaction, the 
magnetic dipole-dipole interaction 17- 24j, which plays 
an important role in domain formation in macroscopic 
samples. This model includes the long-range magnetic 



dipole-dipole interaction between different lattice sites, 
but neglects it within each site[l3|. assuming that it is 
much weaker than the s-wave interaction described. 

It is predicted that the ground state of ^'^Na BEC 
(c2 > 0) is a spin singlet with properties (ni = = 
rt_i = -^/3) [23 contrast with those of mean field pre- 
diction. Soon, Ho and Yip [2^ show that this spin sin- 
glet state is a fragmented condensate with anomalously 
large number fiuctuations and thus has fragile stability. 
The remarkable nature of this fragmentation is that the 
single particle reduced density matrix gives three macro- 
scopic eigenvalues (above) with large number fluctuations 
Ani_o,-i ~ N . In this paper, we consider specifically the 
case of spinor ^^Na condensates and discuss the ground 
state properties and quantum number fluctuations in the 
Zeeman components. 



II. THE MODEL HAMILTONIAN 

The simple two well model Hamiltonian includes H 
Hi + H2 + H12 where 



Hi 



dr 

Pi 



/)2 



Vtrap)'^c 



(1) 



is the in-site spinor BEC Hamiltonian^, the spin- 
dependent 5 interaction is Vs{y) = [a + /3F • F)(5(r) 
with a, (3 characterize the short-range spin-independent 
and spin-changing s-wave collisions, respectively. F — 
al^Yapdp is defined in terms of the 3x3 spin-1 matrices 
Fq^, and a, /3 = 1, 0, — 1 describe the three Zeeman levels 
with repeated indices to be summed over. H2 is identical 
to Hi except for the substitution of subscript 1 by 2 and 
by 

The dipole-dipole Hamiltonian is 



H12 = j dr-^i^lVddi-ri - r2)$/3«'a 



with 



di-d2 3(di • fi2)(d2 • fi2) 



1-12 



1-12 



(2) 



(3) 



2 



fj.Q the magnetic permeability of vacuum, di^2=fl'_FA*_BFi,2 
with gplJ'B the gyromagnetic ratio, and Y12 = ri — r2, 
fi2 — V12/ |ri2|, ri_2 is the coordinate of the 1, 2 site. 

We adopt the single model approximation [H, [g^, HI] 
for each of the two spinor condensates in the nearby sites 
with modes ^'(r) and $(r), i.e., setting 



enhanced by the light-induced optical dipolar interaction 
if one chooses appropriate laser fields to form the poten- 
tial well [2^ [3^, we aim to determine the spin structure 
of the system if the different sites are allowed to inter- 
act with each other through the magnetic dipole-dipole 
interaction. 



1,0,-1 



(4) 



III. THE GROUND STATE PROPERTIES 



with hi [hi) the annihilation operator for the ferromag- 
netic (polar) atoms satisfying [ai,aj] = and ai,aj 

5ij (and the same form of commutations for bi). substi- 
tute into Hi,H2, The spin-independent part can 
be reduced to a constant operator for a fix number of 
atoms iV = iVi + iVo + A^-i , 



A. A brief review of singlet state 

Without the magnetic dipole-dipole interaction or for 
the intra-site pure spin-1 condensate {Hq — CF^), the 
simplest ground state for the F=l spinor ^"^Na conden- 
sates (C > 0) is a spin singlet formed by two spin-1 
particles described as 



(^V^ + C/2 + |$|2)$ = ;,2$ 

here /zi_2 is the mean field energy or the chemical poten- 
tial in the two wells. 

The spin-dependent Hamiltonian finely reduce to 



Hi — CiFi, H2 



C2F2, 



(6) 



with Ci=Pij dr |4'(r)|^ , C2 = /32 / dr |$(r)|^ 

In the one-dimensional double well or optical lattice, 
according to the the vector subtraction, ri2 = ri — r2 
always point at the only one direction. If we choose the 
quantization axis along this direction as z axis, the Vdd 
reduced to 



Vdd = A[Fi-F2-3(Fi-fi2)(F2-fi2)] 

- A[Fi • F2 - 3Fi,F2,] (7) 

with A = i^sImi^.The Hamihonian H12 finely reads 

47r|ri2| -^^ 



H12 = A(Fi • F2 - 3Fi,i^2.) 



with A = A/ dr |^'(r)|^ \<^(r) 
The total Hamiltonian |2(J 



M is 



H 



CiFj 



C2Fi + AFi • F2 - 3AFi,F2, 



(8) 



(9) 



In the absence of long-range magnetic dipole-dipole inter- 
action or external magnetic fields, there is no spin cor- 
relations between sites. For the Rb^^ condensate, the 
ground state in the individual sites favors polarizing all 
the spins to the same direction, therefore they can be 
considered as independent "magnets" whose pseudospin 
vectors point in random directions. But for individual 
Na^'^ condensate, they favors vanishing the total spin in 



\F,m) = 



G|Fi = l,mi) IF2 = l,m2) (10) 



mi ,7712—1,0,-1 

with F = Fi +F2 = is the total spin, m = mi+m2 = 
is the total z component, G is the Clebsch-Gordon coef- 
ficient. 



|F,m) = ^it|0) 



(11) 



The operator = (aj)^ — 2d\a'_^ describe a singlet pair 
creating operator formed by two identical spin-1 bonsons, 
and the ground state of N particles is (A^)^^'^ |0) . G is 
the CG coefficient. The particle density matrix {p)ai3 = 
{dldf}) is 



iV/3 



N/Z 



(12) 



with a, /3=1, 0, — l.This matrix has three equal macro- 
scopic eigenvalues called "superfragmented state" [2^ . A 
weak external magnetic field along z (iJp = CF^ — pFz) 
can break the pairs and polarize the system [3ll |32| with 
the ground state described as 

\F^S,m^S) = (al )'5(it)(Af-S)/2 (^3) 

The particle numbers on the three Zeeman levels are re- 
distributed as 



A^i 



{N + S) {S+\) S 

25-1-3 ^2S^i 
(N-S) (5 + 1) 



No 



N - 



25- 
-5 



25 + 3 

with the 0-component distribution shrink rapidly as S 



each site. As i] = /io(5FMs) /47''|ri2| can be greatly increases. 
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B. A brief review of dimmer state 



For the subspace of exactly one particle per well, the 
Mott-insulator ground state of one-dimensiona l op tical 
lattice has been confirmed to be a dimmer states [40| with 
the form, 



dimer 



56--- 



(14) 



where the state 

*12 = 



= + 1-1,1) -|0,0))i2 

= -^{albl-a\bU~al,bl)\0) 

^^et.|o) 



(15) 



The operator = ajfoj ~ a\b^_^ ~ a^_-^b\ describe a sin- 
glet pair creating operator formed by two spin-1 bonsons 
in the different site. 

For more or at least two particles per well, the ground 
states are more comphcated. 



C. The quaternary-singlet State 

For simplify, we first consider 3AFizF2z = 0, which 
serves as reference case for the complete discussion. 



H = CiFi + C2FI + AFi • F2 
It can be rewrited asl33l-[39l 



H ^aFi + bFi + cF^ 



(17) 



with a = Ci - A/2, 5 = C2 - A/2, and c = A/2,F = 
F1-I-F2 is the total spin operator. The eigenstates of p7)) 
are the common eigenstates for the commuting operators 
Fi,F2,F^, and F^, given by 



\Fi,F2,F,m) 



c'|^;^i;F2,mJ^i>'«i) 1^2,^2) 



(18) 

with the uncoupled basis states |Fi,mi) (|i^2,n^2)) gen- 
erated from equation (jl3p by a repeat using lowering op- 
erator -Fi-, and they can span a Hilbert space of dimen- 
sion [Ni + l){Ni+2)/2^. C is the Clebsch-Gordon 
coefficient. The corresponding eigenenergy is 




FIG. 1: (Color online) The dependence of ground-state order 
parameters on A at fixed values of Ci = 1, C2 = 2, (in the 
unit of |Ci|). Black dashed lines, red short dashed lines, blue 
dot-dashed lines and green solid hnes denote respectively the 
order parameters Fi, F|, Fi ■ F2, and F^. 



Fig. 1 shows the development of the four order param- 
eters 



n 



Fi • F2 

f2 



<Ff >, 

< Fl -^2 >, 

< F^ > 



(20) 



for A (in the unit of |Ci|) with even particles per well ( 
take iVi = 7V2 = 6). We find that in the A < '^'^Cin'^^ 
(16) region, Fj, Fj, and Fi • F2 are all polarized to the max- 
imum with the system being ferromagnetic. In the re- 
gion A g [— Ci — C2,Ci -|- C2], the two sites are essen- 
tially independent for a weak inter-sites dipole-dipole in- 
teraction. This phase is a total spin singlet F^ = de- 
scribed by the direct product of the polar ground state 
(it)A'/2(^t)A^/2 |o) giving rise to F? = F^ = Oand 

Fi • F2 = O.Whcn A > ^gf^^^j'^' they are polarized to 
the maximum but in the opposite directions with F^ = 
and -2Fi • F2 = F? + F| 7^ O.Wc find interestingly that 



in this state the total spin vanishes, while the sites spins 
satisfy Ff = F^ = N{N + 1) . 



The ground state in the region A > 



(2jV-l)C2 

Ci(Ar+i) 



a singlet, with all basis states obeying the condition 
TOi -|- 7712 — 0. All channels of total spin zero have to 
be taken 

into account and we have 



N 



E = aF,{Fi + l) + bF2{F2 + l) + cFiF + l) (19) \N,N,0,0) = ^ C^N"m,;N,-,njN,mi) \N, -m^) 



Given Nj , the allowed values of Fj are Fj — 0,2,4, ■■■ Nj 
if Nj is even; and Fj = 1,3,5, ■■■ Nj ii Nj is odd, satis- 
fying |Fi -F2I i^Fi^Fi+F2. 

We will next consider the special case oi Ni = N2 ^ N 
and for N even. 



(21) 



If we take Ni^N2=2 for example, we find that 



12,2,0,0) = -i=((el2)2-iitst)|o) (22) 



2^/5 
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Compared to the "superfragmented state" {A'^)^^'^ |0) = 
((dj)^ — 2d|o^;^)^/^ |0)[il], this quaternary-singlet states 
has the similar appearance. However it is not the 
direct product of the two conventional singlet pairs, 
|2, 2, 0, 0) ^ it^t |0) . The difference between |iV, N, 0, 0) 
and (it)Af/2(ijt)Ar/2 |o) can be easily find in the phase di- 
agram (Fig.l). 



|26| . they are 



An 



(27) 



For the state ^^/^(gIj)^ |0) , according to the multino- 
mial theorem 



D. The number fluctuations 



2^3 



(28) 



k=0 1=0 



We notice that 



[F?, 612812] ^ 0, [Fi 612812] ^0, 
[F2, 612812] = 



(23) 



and it is easy to understand that F^(6|2)^ |0) = 0. 
In this section, we will talk about three different sin- 
glet states, the quaternary-singlet state |iV, iV, 0, 0) , the 
dimer state (Sj^j)^ |0) and the direct product singlet 
state (i1')^/2(Bl')^/2 |0) , the properties such as 

F^{A+ f' {B+ f^ |0) = 0, F^ \N, N, 0, 0) = (24) 

are shown in Fig.l. 

As a exceptional case, all eigenvectors | A^, mi ) in (fT8|) 
can be expressed in terms of the Fock states [43 , which 
are defined as 



(J) „b)\ _ „0) 



0,±l;i = 1,2 



(25) 



For the state |A^, A^, 0,0) ,we calculate the particle num- 
bers and number fluctuations on the Fock states and flnd 
that the average numbers of atoms in the six components 



are exactly all equal, y/ti y — ytg 
The fluctuations are given explicitly 



A3)\ 



A^/3. 



An, 



0) 



A', 



,0) 



VN^ + 9N 



(26) 



which approximatively satisfy ( An^' 



An^j^l) for large A^f^ 



,0-) 



2 (An 



(J)\ 

; 



as opposed to 2 (Ani) — 

(Ano) = 2 (An_i) for the single species or intra-site sin- 
glet state (Ci,2 > 0) 

The difference between these states are obvious. As 
total spin F vanishes, the number distributions are all 



,0) 



N/3, but the number fluc- 
tuation distribution in these states are quite different, it 
has been shown that for the state Z^/^{A+)'^^ {B+)^^ |0) 



with Cnik = nl/ {ll(k — l)\{n — k)l), we find that the 
state (812)^ |0) can be described by the Fock state 



(1) (1) (1) 

"•1 ,"•() 



(2) ^(2) (2) 



{e\,f |0) = (dJSt - a\bU - aUblYlO) 

N k 

k=0 l=Q 
N k 

= '^^{^if Nl \k - I, N ~ k, I) (E)\l,N-k,k- I) (29) 

k=Q 1=0 

where we have used the property (a^)^ |0) = v^ATjA^). 
We find that the number fluctuations are equally dis- 
tributed, i.e. 



An 



(i) 



An, 



0) 



An 



"0 



N{N + l)/6- N^/9 (30) 



E. The term 3AFi^F2^ 



For the real Hamiltonian the last term plays an 
important role in domain formation and can polarize the 
spin to the same direction. This interaction offers a effect 
extra uniform weak field to the nearby site and breaks 
the singlet states. The ground state can be constructed 
using the quaternary-singlet state |A^, A^, 0, 0) and direct 
product singlet state (it)^/2(Bt)^/2 |0) . 

In the region A e [— Ci — C2, Ci -I- C2], since the spin 
singlet operator commutes with the spin 



[F?,it] 



0, [FlB^] 







(31) 



and it does not change total spin and any spin compo- 
nents but just add two particles. Therefore, we can con- 
struct the unnormalized spin state for N particles [llj: 
first, write down a state with necessary spin for a small 
number of particles; second, apply A'^{B^) as many times 
as needed to get the desired number of particles. We got 

I®) = Z'/\a\fHb\f^iAYHBY' |0) (32) 



5 



with the fixed number in the two sites satisfied Ti 2 = 

(iVi,2 - 5l,2)/2. 

For the region A > if we let |iV,iV,0,0) = 

^i/2(Qt)2Af |o) ^the ground state for the Hamihonian ^ 
is 

|g)=ZV2(at)5(^t)S(Qt)2Ar-2S|o^ (33) 

2iV 1^0 

In the A < CiN ^ region, the system is polarized 
to the ferromagnetic phase, the ground state is 

|P)-Zi/2(at)Wi(5t)W2|o). (34) 



dipole-dipole interaction between nearby wells. We con- 
sider the special case that there are two particles per 
well, and show a new singlet state. In two well model, 
three kinds of spin ground state with total spin vanished 

(F^ = (^(Fi ■¥2)'^'^ = 0) are discussed and can be dis- 
tinguished by the number fluctuations. The final states 
can be constructed by singlet pair creation operator and 
the quaternary-singlet creation operator. 



IV. CONCLUSION 



To summarize, we study the ground spin state of polar 
atoms (^•^Na) in the optical lattice subject to a magnetic 
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